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Abstract. In a quasiclassical framework, we formulate the quadruple ionization by 
single photon absorption of the Coulomb five-body problem. We present the quadruple 
photoionization total cross section of the ground state of Beryllium for energies up to 
620 eV. Our quasiclassical results for energies close to threshold are in agreement with 
the Wannier threshold law for four electron escape. In addition, the agreement of 
our results with a shape formula provides support for the overall shape of our total 
quadruple cross section. Finally, we find that the photon energy where the maximum 
of the total photoionization cross section occurs for single, double, triple and quadruple 
photoionization of H, He, Li and Be, respectively, seems to follow a linear relation with 
the threshold energy for complete break-up of the respective element. 
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1. Introduction 

Multiple ionization of atoms by single photon absorption is a process of great interest 
since it probes the correlated motion of many electrons. Quadruple photoionization 
of the ground state of Beryllium is the most fundamental atomic process involving 
four bound electrons. No other theoretical or experimental study is, to our knowledge, 
currently available regarding escape of four initially bound electrons. 

The ionization processes involving two electron escape can now be treated quite 
accurately and significant progress has been made concerning their understanding [1]. 
Regarding the three electron escape by single photon absorption from the ground state of 
Li, significant advances have been made over the last few years in obtaining total triple 
ionization cross sections both experimentally [2,3] and theoretically [4-6]. However, 
still great challenges remain, particularly concerning differential cross sections [7-10]. 
Triple ionization of the ground state of Li by electron impact, that is, a four electron 
escape process, has been addressed in very recent years. These experimental [11] and 
theoretical [12] studies address very high energies of the impacting electron. 

2. Quasiclassical formulation of quadruple ionization of Beryllium 

Our construction of the initial phase space density p(j) for the quadruple 
photoionization of Beryllium is similar to the double photoionization of He [13] and the 
triple photoionization formulation of Li which has been detailed in [6]. We formulate 
the quadruple photoionization process from the Be ground state (ls 2 2s 2 ) as a two step 
process in the following way: 

a 4+ = a abs P 4+ (1) 

where a abs is the total photoabsorption cross section and P 4+ is the probability for 
quadruple ionization. For a abs we use the experimental data from ref [14]. Our 
formulation accounts for the second step. First, one electron absorbs the photon 
(photo-electron) at time t = t^ s = 0. Through electronic correlations, the energy 
is redistributed, resulting in four electrons escaping to the continuum. We first assume 
that the photon is absorbed by a ls-electron at the nucleus (r x = 0). This latter 
approximation becomes exact in the limit of high photon energy [15]. For Be, the 
cross section for photon absorption from a Is orbital is much larger than from a 2s 
orbital [16]. Hence, we can safely assume that the photo-electron is a Is electron 
which significantly reduces the initial phase space to be sampled. We also neglect 
antisymmetrization of the electrons in the initial state. We denote the photo-electron by 
1, the other Is electron by 2 and the two 2s electrons by 3 and 4, respectively. Following 
photon absorption, we model the initial phase space distribution of the remaining three 
electrons, Is and two 2s, by the Wigner transform [17] of the corresponding initial 
wavefunction ip{v\ = 0, r2, r3, r^, ), where r» are the electron vectors starting at the 
nucleus. We approximate the initial wavefunction as a simple product of hydrogenic 
orbitals 0f'(rj) with effective charges Zj, to facilitate the Wigner transformation. The 
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Zi are chosen so that they reproduce the known ionization potentials Jj, namely for one 
of the two 2s electrons Z 4 = 1.656 (J 4 = 0.343 a.u.), for the other 2s electron Z 3 = 2.314 
(I3 = 0.669 a.u.), and for the Is electron Z 2 = 3.363 (I2 = 5.656 a.u. )(We use atomic 
units throughout the paper if not stated otherwise.) The excess energy, E, is given by 
E = E w — I with E u the photon energy and / = 14.67 a.u. the Be quadruple ionization 
threshold energy. Given the above considerations, the initial phase space density is 
given by 

p( r/ ) = JK6(T 1 )5(E 1 +I 1 -u) J] W^faviWEi + Ii) (2) 

i=2,3,4 

with normalization constant jY . 

We determine the quadruple ionization probability P 4+ formally through 

P 4 + = lim / di> + exp((t - t abs )^f cl )p(r). (3) 

J r abs 

The projector V 4+ indicates that we integrate only over those parts of phase space that 
lead to quadruple ionization. Jzf c i is the classical Liouville operator which is defined by 
the Poisson bracket {H, }, with H the Hamiltonian of the system. In our case H is 
the full Coulomb five-body Hamiltonian. We propagate the electron trajectories using 
the classical equations of motion (Classical Trajectory Monte Carlo method [18,19]). 
Regularized coordinates [20] are used to avoid problems with electron trajectories 
starting at the nucleus. We label as quadruply ionizing those trajectories where the 
energies of all four electrons are positive, E{ > with % = 1,2,3,4, asymptotically in 
time. 

3. Results 

In figure [IJ we present our results for P 4+ . The very small value of the quadruple 
photoionization probability clearly shows the difficulties involved in obtaining theoretical 
results or experimental measurements of the five-body break-up process. Using (CQ) we 
present in figure [2] the results for the quadruple photoabsorption cross section (black 
circles). Since there are no other results currently available we can not compare our 
results on an absolute scale. Regarding the accuracy of our results, the standard relative 
error is oc 1/ WN qua d [21] with N qua d the number of quadruply ionizing events. For 
energies above 10 eV the N qua d we obtain are a few thousands rendering our results very 
accurate while for energies from 4eV up to and 9 eV, N quac i are at least 1000. 

However, a comparison of our results with a shape function applicable to single- 
photon multiple ionization processes [22] provides some indication of how well our results 
reproduce the shape of the total cross section. This shape function reproduces, by 
construction, the correct behavior of the single photon multiple-ionization cross section 
for excess energies close to threshold (Wannier law [23]) and for large excess energies. 
It depends on two parameters, the position Em and height om of the cross section 
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where a is the characteristic exponent in the Wannier threshold law and x is the excess 
energy E scaled by Em- For the break up of a four electron atom this characteristic 
exponent was found to be 3.331 for a plane configuration where the electrons escape 
on the apexes of a square and 3.288 for a three dimensional configuration where the 
electrons escape on the apexes of a tetrahedron [24]. From the above two values for 
a it is the smallest one, a = 3.288, that dominates the threshold behavior and thus 
the value we use in the shape formula in (j4j). Using the shape formula without any 
fitting parameters but substituting for E M = 120 eV and o~ M = 0.0847 b which are 
the location and height of the maximum cross section from our results we obtain the 
solid curve shown in figure [3 Given that there are no fitting parameters the agreement 
between our numerical results and the shape formula is quite good, providing support 
that the shape of the cross section we obtain is the correct one. The agreement becomes 
even better if instead we use the E M and o~m as fitting parameters to fit our results with 
the shape formula. Doing so, we obtain the values E M = 134.5 eV and a M = 0.0834 
b and the agreement of the fitted shape formula (dashed curve in figure [2]) with our 
results is indeed very good. 

Also, for small values of the excess energy we have fitted our results for the 
quadruple cross section with the formula 



with I the fragmentation threshold energy of beryllium and E w the photon energy. For 
energies from 4 eV to 9 eV our results yield a characteristic exponent a of 3.17 which is 
quite close to the analytical value of 3.288. The difference can probably be attributed 
to the fact that our current results go only down to 4 eV. It is a well known fact that as 
the energy increases the value of a decreases [6]. Thus, most probably, if we were able 
to obtain results for energies below 4 eV a would have been closer to the actual value of 
3.288. Currently however the numerical difficulties involved in obtaining results below 
4 eV are prohibitive. 

At this point, a few comments are in order. Our quasiclassical formulation is 
most accurate close to threshold. It can describe the process accurately up to at most 
intermediate excess energies while it can not account for high energies. For high energies 
the process can only be described quantum mechanically. In the case of triple ionization 
from the ground state of Li our simple product initial state has proven to be very good 
in describing not only the total cross section [6] but also double energy differential cross 
sections up to excess energies of 220 eV, see [10], with / = 203.5 eV for the triple escape 
process. How good the choice of a product initial state for the case of beryllium is 
remains to be seen when more theoretical and experimental results become available 
and a comparison with our results can be made for a range of excess energies. 
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Figure 1. Quadruple photoionization probability as a function of excess energy. 
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Figure 2. Quadruple photoionization cross section as a function of excess energy. The 
black circles are our results. For the solid and dashed lines sec the main text. 



A comparison of the maximum of the single photoionization cross section for H, 
6.3 x 10 6 b, with the maxima for double photoionization of He, 8.76 x 10 3 b [25], triple 
photoionization of Li, 7.47 b, and quadruple photoionization of Be, 0.0847 b, clearly 
shows how rare the quadruple photoionization process is. For the double ionization cross 
section of He we use the results presented in ref [25] while for the triple photoionization 
of Li we use the results presented in ref [6] . Another, interesting feature is the relation 
between the ionization threshold energy and the photon energy = I + Em where 
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Figure 3. Photon energy where the maximum of the cross section for single 
photoionization of H, double photoionization of He, triple photoionization of Li, and 
quadruple photoionization of Be occurs as a function of the threshold energy of the 
respective process 

the maximum of the photoionization cross section of the respective process occurs. For 
the single photoionization of H, Em = eV, for the double photoionization of He, 
E M = 21.8 eV [25], while for the triple and quadruple photoionization of Li and Be E M 
is respectively 52.7 eV and 134.4 eV. The values of 52.7 eV and 134.4 eV were found 
by fitting (TjJ to our results for triple [6] and quadruple photoionization, respectively. In 
figure [3], we plot E u as a function of / for the above four processes. Interestingly, we find 
that for these four processes the ratio of E u /I is constant. Whether there is a physical 
reason underlying this constant ratio is an interesting question for future research. 

A future study of interest for the four electron escape in Be would be to investigate 
whether the collision processes the four electrons follow to escape to the continum can 
still be described as a sequence of three-body subsystems as is the case for Li [9]. 
Another interesting aspect to be investigated is the shape of the inter-electronic angular 
distributions for energies close to threshold. In the case, of the three electron escape 
from the ground state of Li we found that in an energy regime where the characteristic 
exponent a = 2.16 [26] is recovered the angular distribution has a surprising T-shape 
distribution. This T-shape is different than what one would expect from the three 
electrons escaping in the apexes of an equilateral triangle for excess energy zero [26]. 
At zero excess energy, the four electrons escape in the apexes of a tetrahedron [24]. 
It remains to be seen whether a surprising pattern will be found for the four electron 
escape as was the case for the three electrons. 

In conclusion, we have presented the first results for the total cross section for 
the quadruple photoionization of Be. With no other theoretical or experimental results 
currently available we hope that the current results will serve as a benchmark calculation 
for future studies of the Coulomb five-body problem. 
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